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Ramsey  spectroscopy  has  become  a  powerful  technique  for  probing  nonequilibrium  dynamics  of 
internal  (pseudospin)  degrees  of  freedom  of  interacting  systems.  In  many  theoretical  treatments,  the  key  to 
understanding  the  dynamics  has  been  to  assume  the  external  (motional)  degrees  of  freedom  are  decoupled 
from  the  pseudospin  degrees  of  freedom.  Determining  the  validity  of  this  approximation — known  as  the 
spin  model  approximation — has  not  been  addressed  in  detail.  Here  we  shed  light  in  this  direction  by 
calculating  Ramsey  dynamics  exactly  for  two  interacting  spin- 1/2  particles  in  a  harmonic  trap.  We  focus  on 
i-wave-interacting  fermions  in  quasi  one-  and  two-dimensional  geometries.  We  find  that  in  one  dimension 
the  spin  model  assumption  works  well  over  a  wide  range  of  experimentally  relevant  conditions,  but  can  fail 
at  time  scales  longer  than  those  set  by  the  mean  interaction  energy.  Surprisingly,  in  two  dimensions  a 
modified  version  of  the  spin  model  is  exact  to  first  order  in  the  interaction  strength.  This  analysis  is 
important  for  a  correct  interpretation  of  Ramsey  spectroscopy  and  has  broad  applications  ranging  from 
precision  measurements  to  quantum  information  and  to  fundamental  probes  of  many-body  systems. 
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Ramsey  spectroscopy,  a  technique  initially  designed  to 
interrogate  microwave  atomic  clocks,  has  become  an 
important  modern  tool  for  probing  dynamics  of  interacting 
many-body  systems  with  internal  (pseudospin)  degrees  of 
freedom.  Ramsey  specttoscopy  applies  [see  Eig.  1(a)]  two 
strong  resonant  pulses  to  a  system  initially  prepared  in  a 
well-defined  pseudospin  state,  separated  by  a  dark  time  of 
free  evolution.  The  first  pulse  initializes  the  pseudospin 
dynamics  by  preparing  the  system  in  a  nontrivial  super¬ 
position  of  eigenstates;  i.e.,  it  introduces  a  quantum  quench 
[1].  The  second  pulse  reads  the  coherence  or  correlations 
developed  during  the  dark  time.  Recently,  Ramsey  spec¬ 
troscopy  has  been  proposed  for  extracting  real-space  and 
time  correlations  [2-6],  characterizing  topological  order 
[7,8],  measuring  spin  diffusion  dynamics  in  bosonic  [9-14] 
and  fermionic  systems  [15-18],  and  as  a  means  to  probe 
many-body  interactions  in  atomic,  molecular,  and  trapped 
ion  systems  [19-27]. 

Generally  speaking,  Ramsey  spectroscopy  measures  the 
collective  pseudospin  and  traces  out  other  external  degrees 
of  freedom  involved  during  the  free  evolution.  In  most 
atomic  setups  the  latter  are  associated  with  motional 
degrees  of  freedom  in  the  harmonic  trapping  potential 
and/or  lattice  potential  confining  the  atoms.  The  external 
degrees  of  freedom  can  affect  the  spin  dynamics  in  a 
nontrivial  way,  however.  A  great  simplification  could  be 
gained  if  it  were  possible  to  decouple  the  motional  and  spin 
degrees  of  freedom,  and  reduce  the  many-body  dynamics 
down  to  those  extracted  from  a  pure  interacting  spin 
model.  Evidence  that  this  scenario  is  possible,  even  far 
from  quantum  degeneracy,  has  been  reported  in  recent 


experiments  [9,13,14,25-28],  where  the  observed  spin 
dynamics  corresponded  to  those  of  a  pure  spin 
Hamiltonian.  These  observations  are  opening  a  path  for 
the  investigation  of  quantum  magnetism  in  atomic  systems 
without  the  need  for  ultralow  temperatures.  It  is  thus 
important  to  determine  the  parameter  regime  in  which  a 
pure  interacting-spins  picture  is  valid. 

In  this  Eetter  we  provide  insight  on  the  validity  of  a  pure 
spin  model  description  of  Ramsey  spectroscopy  by  perform¬ 
ing  exact  calculations  for  fermions  with  5-wave  interactions 
and  an  internal  pseudospin- 1/2  degree  of  freedom,  confined 
in  quasi- ID  and  quasi-2D  harmonic  traps.  We  show  that  the 
large  degeneracy  of  the  harmonic  oscillator  spectrum  can 
limit  the  validity  of  the  spin  model  to  time  scales  less  than  the 
inverse  interaction  strength,  due  to  resonant  collisionally 
induced  excitation  of  spatial  modes  [see  Eig.  1(b)].  Cold 
atom  experiments  are  protected  from  this  problem  if  the 
temperature  is  high  enough  that  atoms  probe  the  actual 
Gaussian  shape  of  the  potential  which  breaks  the  harmonic 
spectrum  degeneracy.  This  was  shown  to  be  the  case  for 


FIG.  1  (color  online),  (a)  Ramsey  spectroscopy  of  two  inter¬ 
acting  spin- 1/2  particles,  (b)  In  a  harmonic  trap  the  spectmm 
degeneracy  allows  near-resonant  mode-changing  collisions 
coupled  to  the  spin  dynamics. 
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example  in  Refs.  [9,25,26]  where  a  pure  spin  model  well 
described  the  experimental  observations.  At  very  low 
temperatures,  Pauli  blocking  can  also  prevent  mode  chang¬ 
ing  collisions,  as  recently  observed  in  Ref.  [29].  However, 
the  degeneracy  is  a  concern  for  intermediate  temperatures  at 
which  the  set  of  populated  levels  is  effectively  harmonic. 
Here  we  show  that  surprisingly,  in  two  dimensions  and  to 
first  order  in  the  interaction  strength,  the  full  two-particle 
dynamics  can  be  described  in  terms  of  an  effective  spin 
model  with  appropriate  parameters.  Our  two-body  calcu¬ 
lations  are  not  only  a  first  step  towards  understanding  the 
interplay  between  spin  and  particle  motion  in  generic  many- 
body  ensembles,  but  are  also  directly  applicable  to  optical 
clocks  that  interrogate  an  array  of  ID  tube-shaped  traps, 
each  with  fewer  than  three  atoms  [28,30,31]. 

Physical  situation. — Consider  two  fermions  with  internal 
degrees  of  freedom  {t>  ^l'}  corresponding,  for  instance,  to 
the  —  ^Pq  electronic  levels  in  alkaline-earth-based 
optical  lattice  clocks,  and  assume  their  interactions  are 
primarily  described  by  an  5-wave  pseudopotential.  The 
atoms  are  also  illuminated  by  a  laser  beam  detuned  by  5  = 

Wl  —  coq  from  the  atomic  transition  qjq,  with  wave  vector  k 
and  bare  Rabi  frequency  Q.  The  two-particle  Hamiltonian 
is  then  given  by  H{xi,X2)  =  E/=i,2^Lft)  +  ^d(^h^2): 

+H.C., 

Hd{xi,X2)  =  H,p{xi)  +  H ,p{x2)  +  gPsS{r)  -^r.  (1) 

Here,  H\^{xi)  describes  the  atom-laser  interaction:  o-t  is  the 
spin  raising  operator  acting  on  atom  i,  and  H.c.  is  the 
Hermitian  conjugate.  Hsp{xi)  =  — /i^/(2M)V?  -|-  V(3c,)  -|- 
{ha)(^l2)a]  is  the  single  particle  Hamiltonian  with  an 
external  potential  V,  assumed  for  simphcity  to  be  inde¬ 
pendent  of  the  internal  state  and  separable.  H^p{xi)  has 
eigenfunctions  <Pn(xi)  and  eigenenergies  £„  with 
n  =  {wj.,  Wy,  Wj}.  M  is  the  particle’s  mass  and  the 
Pauli  matrix,  r  =  3ci  —  3c2  is  the  relative  coordinate, 
g  =  Anh^al^ /M  and  nj'*'  the  3D  5-wave  scattering  length. 
Pj  =  |5)(5|  is  the  projector  into  the  singlet  state, 
I5)  =  (l/\/2)(|t|)  —  lit))-  Only  fermions  in  the  singlet 
state  can  interact,  while  spin  triplet  states,  =  Ij.],), 

kit)  =  (W2)(|tl)  +  lit)),  and  =  Iff)  cannot 
experience  5-wave  interactions. 

The  spin  model. — The  assumptions  of  the  spin  model 
are:  if  there  are  no  degeneracies  in  the  two-atom  non¬ 
interacting  spectrum,  i.e.,  [E^  -!-£'„)  =  (£„,/  -|-  E„i) 

occurs  only  for  (m,n)  =  (m',n')  or  (m,n)  =  (n',m'), 
and  interactions  are  treated  as  a  perturbation,  scattering 
processes  that  change  the  single-particle  modes  become  off 
resonant  and  atoms  remain  frozen  during  the  dynamics.  In 
this  case  interactions  are  diagonal  in  the  single -particle 


basis  and  for  particles  in  modes  (m,  n)  they  are  fully 
characterized  by  the  interaction  energy 

hU'\'l  =  g  J  (2) 

Eermions  with  5- wave  interactions  in  one  dimension. — 
We  begin  with  the  case  of  two  atoms  tightly  confined 
transversally  in  their  ground  state  and  with  dynamics  only 
along  the  z  direction,  where  they  experience  a  ID  harmonic 
trapping  potential  with  angular  trapping  frequency  tw^. 
The  two  atoms  are  initially  prepared  in  the  state 

(1/x/2)(|wi,W2)  -  l«2.«i))ku)- 

The  atoms  are  assumed  to  be  in  the  Eamb-Dicke  regime, 
with  Eamb-Dicke  parameter  rj  =  k^afigfs/2  ^  1.  a^,,  = 
is  the  harmonic  oscillator  length,  and  k^  the 
projection  of  the  probe  laser  wave  vector  along  z.  Mode 
changes  during  the  laser  interrogation  can  be  suppressed 
if  the  laser  detuning  from  the  atomic  transition  3  and  the 
bare  Rabi  frequency  Q  satisfy  S,  tjQ.  ^a>^.  In  this  regime 
the  mode  dependence  of  the  Rabi  frequencies  is  = 
Qe-('7V2)£0('^2^  [32].  The  Hamiltonian  in  the  rotating 
frame  of  the  laser  [21-23,25]  under  the  spin  model 
approximation  can  be  written  as  H'shf^  = 
where 

L  2 

Ho"”"  =  2hu':^'^"^P,  -  hSs^.  (3) 

acts  only  during  the  two  laser  pulses,  and  acts 
only  during  the  dark  time.  Here,  5;,.^ ^  -|-  02^’^) {2 

are  collective  spin  operators  and  +  Q„^)/2  is 

the  mean  Rabi  frequency.  AQ''>  ''2  =  —  Q„^)/2  arises 

from  the  excitation  inhomogeneity  and  can  transfer  some  of 
the  initial  triplet  population  to  the  singlet,  allowing  inter¬ 
actions.  The  interaction  energy  =  hU'"’^  in  Eq.  (2) 

with  n  =  {0, 0,  Wl}  and  m  =  {0, 0,  W2}.  We  can  ignore  the 
detuning  and  interactions  during  the  laser  pulses  if  the 
pulses  are  short  compared  to  the  time  scales  set  by  those 
energies.  We  also  ignore  single-particle  energies  which  are 
constants  and  do  not  contribute  to  the  dynamics. 

The  spin  model  assumptions  break  down  in  a  harmonic  trap 
due  to  the  degeneracy  of  the  noninteracting  two-atom  spectmm: 
even  weak  interactions  can  transfer  atoms  initially  in  modes 
{ni,n2}  to  the  various  degenerate  configurations  {n\ + 
k,  W2  —  k}  (for  integer  k)  during  the  dynamics.  To  account 
for  these  mode  changes,  we  take  advantage  of  the  exact 
eigenfunctions  and  eigenvalues  of  Hd{xi,X2)  in  Eq.  (1)  for 
two  atoms  with  5 -wave  interactions  in  a  harmonic  trap  [33]. 
These  solutions  exploit  the  separability  of  the  Hamiltonian  in 
the  center-of-mass  coordinate  R  and  relative  coordinate  r. 
There  is  no  degeneracy  in  the  relative  coordinate  degree  of 
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freedom.  See  Supplemental  Material  [34]  for  straightforward 
expressions  for  the  change  of  basis.  Equivalent  expressions  are 
given  in  [35]. 

Ramsey  dynamics  in  the  spin  model  approximation. — 
Denoting  r  the  Ramsey  dark  time,  the  population  difference 
between  the  two  spin  states  measured  after  the  second  pulse 
takes  the  generic  form 


{s^){'t)  =  A(t)  cos((5t)  +  B{t)  sin((5T)  +  C(t).  (4) 

A(t),  B{t),  and  C(t)  have  the  form  A(t)  =  +  /2, 

B{t)  =  /2(t)/3,  C(t)  =  /3(t)/4  +/5.  where  /,(t)  depend 
on  the  dark  time  physics,  and  /,■  are  independent  of  the 
dark  time  physics  and  depend  only  on  the  laser  pulse 
quantities  (see  Supplemental  Material 

[34]).  and  0"*’"^  =  with  fi,2 

the  pulse  durations.  In  the  spin  model  approximation,  the 
dark  time  functions  depend  simply  on  interactions: 

;sm  ^  pm  ^  cOs{u'l'f^T),  =  sin(M"‘j^^'’^T). 

Ramsey  dynamics  in  the  weakly  interacting  regime 
m^). — ^Eor  weakly  interacting  atoms  (Mj.j’  ^  m,), 
we  are  able  to  write  the  dynamics  (beyond  the  spin  model 
approximation)  in  a  closed  analytic  form  [34].  These  expres¬ 
sions  for  the  dynamics  are  exact  for  times  t 


11 1 


/exact  ^  /exact  =  2  ^  \dn'; 


'i2|2 


COS 


n,.=0.even 


'AE,{n,) 


h 


11  { -\-ll2 


/exacts  2  ^ 

i2r=0,even 


sm 


h 


(5) 


Here,  are  the  change  of  basis  coefficients  defined 

in  Supplemental  Material  [34].  Comparing  Eq.  (5)  to  the  spin 
model  solution,  we  see  the  single  frequency  in  the  spin- 
model  dynamics  gets  replaced  by  a  sum  over  many  frequencies 
AEgipi^jh  in  the  exact  dynamics.  These  frequencies  are 
associated  with  the  first  order  correction  of  the  eigenenergies 
due  to  interactions  [33]:  AE^(n^)  =  hu^.^^{Y{n,./2  +  1/2)/ 
y^r(wA2  +  1))(1  +  C>(m|’|/cOj)).  The  many  frequencies 
that  appear  come  ifom  the  resonant  mode-changing  processes. 
States  with  odd  n^  do  not  experience  ^-wave  interactions  and 
do  not  contribute. 

When  we  compare  the  exact  dynamics  to  those  predicted  by 
the  spin  model  we  find  that  they  agree  for  short  times, 
1.  The  spin  model  fails  at  longer  times,  however, 
when  leakage  of  the  population  to  other  modes  in  the 
individual-particle  coordinate  basis  becomes  significant  [See 
Eig.  2(a)].  This  is  reflected  in  the  behavior  of  the  angular 
frequency  shift  Aa»(T) — an  important  quantity  for  atomic  clock 
experiments — defined  as  A(u(t)t  =  — arctan[fi(T)/A(T)], 
which  is  the  observed  change  in  the  atomic  transition  due  to 
interactions  [see  Eig.  3(a)].  The  failure  of  the  spin  model  at 


0  1000  2000  3000  4000  5000 

time  (1/fcjj.) 


FIG.  2  (color  online).  Ramsey  dynamics  [seeEq.  (4)]  withd  =  0: 
(a)  ID  spin  model,  exact  solution,  and  projection  of  population  onto 
initial  mode  (here  n  1  =  10and«2  =  0),  with  Mil*  «  0.2(Uj..  Dephas¬ 
ing  of  the  exact  dynamics  results  from  mode  changes,  (b)  Thermal 
averages  in  two  dimensions:  spin  model  vs  effective  spin  model,  at 
different  temperatures,  with  mIJ*  » 0.04(Wj^.  For  both  figures: 
01  =  02  =  7i/3,  with  thermally  averaged  inhomogeneity 
( Afl) / (fl)  =  0.3.  0i  =  Q.ti  are  bare  pulse  areas. 

times  longer  than  the  inverse  interaction  strength  limits  its 
applicability  to  model  the  new  generation  of  atomic  clocks  that 
use  ultra  coherent  lasers  [36,37],  allowing  interrogation  times 
exceeding  a  few  seconds.  A  spin  model  treatment  will  be 
insufficient  when  conditions  are  such  that  the  atoms  see  an 
almost  purely  harmonic  potential. 

Ramsey  dynamics  in  the  strongly  interacting  regime 
(mJ|.'J'  >  oij). — The  spin  model  fails  when  mII*  >  to^.  To 
maintain  the  separation  between  interaction-induced  effects 
and  laser-induced  effects,  we  imagine  interactions  set  to  be 
weak  during  the  laser  pulses  and  suddenly  increased  after 
the  first  pulse  using,  for  example,  a  Eeshbach  resonance 
[25,38,39]  [40].  For  this  situation,  we  can  solve  for  the 
dynamics,  given  an  initial  pair  of  modes,  although  there  is 
no  closed  form  solution  (the  dark  time  functions  /,  (t)  are 
more  complicated,  but  the  laser  dependence  through  fj 
remains  the  same  as  in  the  previous  cases).  We  find  that,  in 
the  limit  of  strong  interactions  (mIJ*  »  co^),  the  population 
imbalance  exhibits  periodic  oscillations  at  the  axial  trap¬ 
ping  frequency  cOj,  in  contrast  with  the  spin  model 
prediction  of  much  faster  oscillations  at  the  interaction 
frequency  (see  Supplemental  Material  [34]).  The  frequency 
shift  (proportional  to  this  oscillation  frequency),  saturates 
to  a  value  on  the  order  of  co^,  instead  of  increasing  without 
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FIG.  3  (color  online).  Thermally  averaged  frequency  shifts: 
(a)  ID  spin  model  and  exact  solution  vs  population  excitation 
fraction  (number  of  atoms  in  f  divided  by  the  total  number  of 
atoms)  after  the  first  pulse,  at  intermediate  and  long  times.  Here, 
«  0.2aij.  (b)  Frequency  shifts  for  2D  spin  model  vs 
effective  spin  model,  with  «  O.Odcuj^.  For  both  figures: 

T  =  lOShw^/kg,  9  =  0\  =  02,  (Wj  =  cux  =  700  x  In  Hz,  and 
thermally  averaged  inhomogeneity  (AD)/(fl)  =  0.3.  9i  =  Qti 
are  bare  pulse  areas. 

bound.  These  results  reflect  the  fact  that  for  strong 
interactions  (unitarity),  the  fermions  maximally  repel  each 
other,  and  the  trap  energy  becomes  the  only  relevant  energy 
scale  in  the  system.  This  behavior,  expected  to  be  a 
universal  result,  should  apply  even  in  the  many-body  case 
as  seen  in  Refs.  [19,20]. 

References  [21,25]  showed  that  ^-wave  frequency  shifts 
can  be  cancelled  by  setting  the  second  pulse  area  to 
9"''"^  =  11  j 2.  This  result,  obtained  using  the  spin  model, 
survives  the  inclusion  of  resonant  mode  changes  even  for 
strong  interactions  during  the  dark  time,  since  the  depend¬ 
ence  of  the  dynamics  on  the  functions  /,,  and  thus 
remains  the  same  even  when  interactions  are  strong. 

Fermions  with  s-wave  interactions  in  two  dimensions. — 
For  an  anisotropic  2D  harmonic  potential  with  no  acci¬ 
dental  degeneracies,  the  treatment  will  be  similar  to  the  ID 
case.  An  isotropic  2D  harmonic  potential,  however,  is  more 
difficult  to  treat,  due  to  the  large  degeneracy.  In  two 
dimensions  the  spin  model  remains  the  same  as  Eq.  (3), 
with  populated  modes  now  w,-  =  {n^i,nyj,0},  and  inter¬ 
action  energy  in  Eq.  (2).  To  go  beyond  the 


spin  model  we  use  polar  relative  coordinates  to  eliminate 
much  of  the  degeneracy.  For  noninteracting  particles,  the 
eigenfunctions  can  be  parametrized  by  quantum  numbers  n 
and  m,  with  energy  E  =  hco_\_{2n  +  \m\  +  1)  and  angular 
momentum  component  =  hm,  where  is  the  2D 
oscillator  frequency.  S-wave  interactions  only  affect 
states  with  m  =  0,  and  this  subset  of  states  contains  no 
degeneracy  (other  than  the  center-of-mass  degeneracy). 
To  first  order  in  perturhation  theory  the  interaction  energy 
shift  is  independent  of  the  radial  quantum  number  n\ 
AE  =  {g/2s/27i^/^a^a\){l  +  0{u]^y co _[_)),  where  and 
Qx  are  the  oscillator  lengths  along  the  tightly  confined  z 
direction  and  the  weakly  confined  x  and  y  directions, 

respectively,  and  in  2D  with  wj  =  (1,0,0), 

W2  =  (0,0,0).  This  result  is  striking:  despite  the  large 
degeneracy  in  two  dimensions,  each  interacting  state  with 
m  =  0  receives  the  same  energy  shift  to  first  order  in 
perturhation  theory,  and  accumulates  the  same  phase  during 
the  dark  time.  An  effective  spin  model,  with  diagonal 

matrix  element  2hu'^'^^  replaced  hy  AE,  will  he  exact  for 
the  m  =  0  states,  to  first  order  in  the  interaction  strength. 
We  can  replace  in  Eq.  (3)  with 

=  AEE,„=o  -  fiSK,  (6) 

where  .P,„=o  projects  onto  interacting  states  with  m  =  0. 
For  a  properly  symmetrized  initial  state  ^'ni.n2(^i  Az)  in 
modes  (nj,  n2),  we  denote  the  fraction  of  the  population 
with  m  =  0  in  the  relative  coordinate  hy  which 

can  he  calculated  as  =  j 

8{x\  —X2)  =  (4v/2;r^/^/p)hM'||”^,  where  is  the  2D 
interaction  energy  calculated  from  Eq.  (2).  The  dark  time 
dynamics  of  this  effective  spin  model  are  simple: 
/f  ”  =  /f”  =  (  1  -  P'"=\,.n2)  +  ^n.Tn”  COs(AET//t),  /f  ”  = 
P'"=S,  sin(AET//i). 

In  the  original  spin  model,  is  used  as  the 

interaction  energy.  We  see  that  this  parameter  appears  in 
the  effective  spin  model  to  quantify  the  population  of 
interacting  modes  (En,7n2)’  instead  of  their  energy.  This 
dramatic  result  is  seen  in  Fig.  2(h),  comparing  thermal 
averages  of  the  previously  implemented  spin  model  with 
the  new  effective  spin  model.  Oscillations  during  the 
dynamics  remain  at  the  same  frequency  AE  at  higher 
temperatures,  hut  the  amplitude  of  the  oscillations,  propor¬ 
tional  to  decreases.  The  previously  implemented 

spin  model,  on  the  other  hand,  predicts  smaller  interaction 
energies  (slower  oscillations)  at  higher  temperatures.  The 
frequency  shift  predicted  hy  the  original  spin  model  is  only 
valid  at  short  times  [see  Fig.  3(h)]. 

Summary  and  outlook. — ^We  test  the  validity  of  a  spin 
model  treatment  for  Ramsey  spectroscopy  with  exact 
calculations  for  two  pseudospin- 1/2  fermions  in  a  harmonic 
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trap.  In  one  dimension  the  spin  model  treatment  breaks 
down  for  dark  times  on  the  order  of  the  inverse  interaction 
strength,  and  for  strong  interactions.  In  two  dimensions 
we  find  an  effective  spin  model  which  is  exact  to  first  order 
in  perturbation  theory,  and  whose  dynamics  can  be  quite 
different  from  those  predicted  by  a  spin  model  treatment. 
Euture  theoretical  treatments  of  interacting  systems  probed 
by  Ramsey  spectroscopy  must  take  these  effects  into 
account  to  correctly  describe  dynamics  outside  of  the 
short-time  and  weakly  interacting  regimes. 
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